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The phenomenological Maxwell field is quantized for arbitrarily space- 
and frequency-dependent complex permittivity. The formalism takes account 
Q"^ ■ of the Kramers-Kronig relation and the dissipation-fluctuation theorem and 

On . yields the fundamental equal-time commutation relations of QED. Applica- 

tions to the quantum-state transformation at absorbing and amplifying four- 
port devices and to the spontaneous decay of an excited atom in the presence 
of absorbing dielectric bodies are discussed. 
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O \ Quantization of the electromagnetic field in dispersive and absorbing dielectrics requires 

a concept which is consistent with both the principle of causality and the dissipation- 

0\ • fluctuation theorem and which necessarily yields the fundamental equal-time commutation 

relations of QED. In order to achieve this goal, several approaches are possible. The micro- 
scopic approach starts from the exact Hamiltonian of the coupled radiation-matter system 
and integrates out, in some approximation, the matter degrees of freedom to obtain an ef- 
C ' fective theory for the electromagnetic field. Since the procedure can hardly be performed for 

arbitrary media, simplified model systems are considered. A typical example is the use of 
O^. harmonic-oscillator models for the matter polarization and the reservoir variables together 

with the assumption of bilinear couplings [|[] . In the macroscopic approach, the phenomeno- 
logical Maxwell theory, in which the effect of the medium is described in terms of constitutive 






x 



Vh ' equations, is quantized. Since this concept does not use any microscopic description of the 

medium, it has the benefit of being universally valid, at least as long as the medium can be 
regarded as a continuum. 

Here we study the problem of quantization of the phenomenological Maxwell theory 
for nonmagnetic but otherwise arbitrary linear media at rest, starting from the classical 
Green function integral representation of the electromagnetic field. The method was first 
established for one-dimensional systems [Q and simple three-dimensional systems |]^ and 
later generalized to arbitrary inhomogeneous dielectrics described in terms of a spatially 
varying permittivity which is a complex function of frequency [Q. 

In Sec. H we briefly review the quantization scheme and give an extension to anisotropic 
dielectrics (including amplifying media), which complete the class of nonmagnetic (local) 
media. In Sec. |T| we apply the method to the problem of quantum-state transformation 



at absorbing and amplifying four-port devices, and in Sec. IV we give an application to the 



problem of spontaneous decay of an excited atom in the presence of absorbing bodies. 



II. QUANTIZATION SCHEME 

Let us first consider the electromagnetic field in isotropic dielectrics without external 
sources. The (operator- valued) phenomenological Maxwell equations in the temporal Fourier 
space read 

V-B(r,cj) = 0, V xE{r,uj) = iujB{r,uj), (1) 

V ■ [eoe(r, a;)E(r, u)] = p(r, uj), V x B(r, u) = ~i{uj/c^)e{r, uj)E{r, u) + /ioj(r, uj). (2) 

From the principle of causality it follows that the complex-valued permittivity e(r, u) = 
eR{r,uj) + iej{r,uj) satisfies the Kramers-Kronig relations. Hence, it is a holomorphic func- 
tion in the upper complex frequency plane without poles and zeros and approaches unity in 
the high-frequency limit. Consistency with the dissipation-fluctuation theorem requires the 
introduction of an operator noise charge density p(r, u) and an operator noise current den- 
sity j(r, cj) satisfying the equation of continuity. Quantization is performed by introducing 
bosonic vector fields ffr, cj), 



j(r, uj) = uj^heoeiir, uj)/tx f (r, uj), (3) 

which play the role of the fundamental variables of the theory. All relevant operators of 
the system such as the electric and magnetic fields and the matter polarization can be 
constructed in terms of them. For example, the operator of the electric field is given by the 
integral representation 

^fc(r) = i/ioy/ieo/TT / duj j d^r uj^^Jej{r',uj)Gkk'{r,r',uj)fk'{r',uj) + H.c, (4) 

with Gfcfc'l'') '''5'-^) being the classical dyadic Green function. This representation together 
with the fundamental relation 

/ d^s{uj/cYej{s,uj)G^k{r,s,uj)G*f^{r',s,uj) =lmGij{v,r',uj), (5) 

which follows directly from the partial differential equation for the dyadic Green function, 
leads to the equal-time commutation relation [^ 

eoEk{r),Bi{r') = (h/n) eimk'd^, du {uj/c^)Gkk'{r,r',u). (6) 

Using general properties of the Green function, it can be shown @ that Eq. (^) reduces, for 
arbitrary e{r,uj), to the well-known QED commutation relation 

eo4(r), Bi{r')] = iheuimdt5{v - r') (7) 

The extension to anisotropic and amplifying media is straightforward, since we may 
assume the medium to be reciprocal, so that the permittivity tensor ejj(r,co') is necessarily 
symmetric. In particular, ejj(r, tu) can be diagonalized by an orthogonal matrix Ofc/(r, cj). 
With regard to amplifying media, we note that amplification requires the role of the noise 
creation and annihilation operators to be exchanged. The calculation then shows that the 
fundamental relation (Rl) can be generalized to 



j.(r,^) = u^heo/n [-fij{r,uj)fj{r,u)+-f+j{r,u)f]{r,u)\ , (8) 

^^ 7,^(r,^) = Oikir,uj)^J\iMiir,uj)\ Oij\r,uj) [±eH/(r,^)] , (9) 

eiji{r,uj) = 6ijef\r,uj) = 0:ii}{r,uj)eM i{r,uj)Oij{r,uj). (10) 

Equation (H) completes the quantization scheme for the electromagnetic field in arbitrary 
linear, nonmagnetic (local) media. 



III. QUANTUM-STATE TRANSFORMATIONS BY ABSORBING AND 
AMPLIFYING FOUR-PORT DEVICES 

Let us first apply the theory to the problem of quantum-state transformation at absorbing 
and amplifying four-port devices such as beam-splitter-like devices. Specifying the formulas 
to the one-dimensional case for simplicity and rewriting the integral representation (0) in 
terms of amplitude operators %(a;) and &j(a;) for the incoming and outgoing waves (j = 1, 2), 
the action of an absorbing device can be given by the (vector) operator transformation 

h{uj) = T{uj)k{uj) + A{uj)giiu), (11) 

where gj{uj) are the operators of device excitations and T{uj) and A{uj) are the characteristic 
transformation and absorption matrices of the device given in terms of its complex refractive- 
index profile 1^. Note that aj{u) and gj{uj) are independent bosonic operators. Further, 
it can be shown that the relation T{uj)T^{uj) + A(ti;)A+(ti;) = I is satisfied, which ensures 
bosonic commutation relations for hj{uj). In order to construct the unitary transformation, 
we introduce some auxiliary (bosonic) device variables hjiu), combine the two- vectors a(ti;) 
and g(a;) to the four- vector a(a;), and accordingly b(ti;) and h.{uj) to ^{uj). The four- vectors 
cx.{ijj) and I3{oo) are related to each other as 

^(cu) = A(cj)a(cj), A(cu) G SU(4). (12) 



Introducing the positive Hermitian matrices C(c<j) = wT(ci;)T+(a;) and S{uj) = JA{uj)A+{uj)^ 
the four- matrix A (a;) can be written in the form 

(A = 1). The input-output relation (0) can then be expressed in terms of a unitary operator 
transformation /3(a;) = JJ'^6l{uj)JJ . Equivalently, U can be applied to the density operator of 
the input quantum state pin[a(ct;), q;'(c<j)], and tracing over the device variables yields 

^(F.eld) ^ ^^(Dcvice) {p^„ [a+(^) A(^), A^(^)«t (^)] } . (14) 

To give an example, let us consider the case when one input channel is prepared in an 
n-photon Fock state and the device and the second input channel are left in vacuum, i.e.. 
Pin = |ra, 0, 0, 0)(?7,, 0, 0, 0|. Applying Eq. (|T3|), after some algebra we derive for the density 
operator of the z-th output channel 

«r' = £(';)m.i^'(i-miiT"i*)w. (15) 

Next, let us assume that the two input channels are prepared in single-photon Fock states, 
i.e.. Pin = |1, 1, 0, 0) (1, 1, 0, 0|. We derive for the density operator of the z-th output channel 

pfui? = [1 - \T.i\\l - \T.2?) - \T.2\\l - \Ta?)\ |0)(0| 

+ (|T,i|' + |T,2p - 4|T,inT,2|') |1)(1| + 2|T,inT,2n2)(2|. (16) 

The extension to amplifying devices is straightforward. One has to replace the annihila- 
tion operators Qjioj) in Eq. ( |TT| ) by the corresponding creation operators gj{uj)- This leads 
again to an input-output relation of the form (|12D but with A = — 1 in Eq. ([T3|), the matrix 
A{uj) being now an element of the noncompact group SU(2,2). 



IV. SPONTANEOUS DECAY NEAR DIELECTRIC BODIES 

Spontaneous decay of an excited atom is a process that is directly related to the quantum 
vacuum noise, which in the presence of absorbing bodies is drastically changed and so is the 
rate of spontaneous decay, because of the additional noise introduced by absorption. To 
study a radiating (two-level) atom in the presence of dielectric media, we start from the 
following Hamiltonian in dipole and rotating wave approximations: 

dh dujhujf{r,uj)-i{r,uj) + yhuJo,A^^- iu;2i^2iA(+)(rA) • dsi + H.c. . (17) 

Here, the atomic operators A^a' = |tt)(tt'| are introduced, and A*^+)(rA) is the (positive- 
frequency part of the) vector potential (in Weyl gauge) at the position of the atom. Note 
that the first term in Eq. ( [T7| ) is the (diagonal) Hamiltonian of the system that consists of 
the electromagnetic field and the medium (including the dissipative system) and is expressed 
in terms of the fundamental variables f(r,u;). Solving the resulting equations of motion in 
Markov approximation, the well-known Bloch equations for the atom are recognized, where 
the decay rate is given by [0] 

r = 2u\fikfik' / (heoc^) lmGkk'{rA,rA,ujA) (18) 

[/ifc = {d2i)k, ^A = '^2i]- Note that from Eq. (^) together with Eq. (|^) it follows that 

lmGkkir,r\u)6{uJ-uj') = 7Teocy{hu'){0\[E,{r,u;),E!k'{r',u;')]\0) (19) 

in full agreement with the dissipation-fluctuation theorem. 

Equation ([l8| ) is valid for any absorbing dielectric body. For example, when the atom is 
sufficiently near to an absorbing planar interface, then purely nonradiative decay is observed, 
with § 




\e{iUA) + l\H2LUAz)^ 



where z is the distance between the atom and the interface, and Fq is the spontaneous 
emission rate in free space (for a guest atom embedded in an absorbing dielectric, see |^). 
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